Momentum-dependent light scattering in a 2D Heisenberg antiferromagnet 
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Motivated by the achievements of the X-ray scattering technique, we analyzed the profile of 
the light scattering intensity R{q, uj) at a finite g in a 2D Heisenberg antiferromagnet. Previous 
Raman scattering studies at g = identified the two-magnon peak in Big scattering geometry. We 
found that the Big peak disperses downwards at a finite q, and its intensity increases, reaching its 
maximum at g = go = (0, n) and symmetry related points. In addition, the intensity in the Aig 
geometry becomes non-zero at a finite q, and also displays a two-magnon peak which gains strength 
and disperses to larger frequencies with increasing q, and reaches its highest intensity at go- We 
found that the profile of R{qo,uj) is equivalent in Aig and Big geometries. 

PACS numbers: 78.70.Ck, 72.10.Di, 78.30.-j, 75.30.Ds 



Raman scattering is a powerful probe of magnetic and 
electronic correlations in interacting electron systems. 
Raman studies of antiferromagnetic parent compounds 
of the cuprates in early 90's provided the first estimates 
of the Heisenberg exchange integrals in La2CuOi and 
y^zCsOe- The values of J - 120 - UOmeV were ex- 
tracted from the positions of the Big two-magnon peak 
which, according to the theory, are located at approxi- 
mately 2. 7 J [ij. The values of J extracted from the Ra- 
man studies were later confirmed by neutron scattering 
measurements 1^. 

The traditional framework for the understanding of the 
two-magnon Raman scattering in antiferromagnets has 
been the semi-phenomenological Loudon-Fleury model 
for the interaction of light with spin degrees of free- 
dom U . The model assumes the existence of the matrix 
element Mr for the direct process in which incident pho- 
ton with energy Ui comes in, an outgoing photon with 
energy cuf comes out, and two magnons are excited with 
energies f2(fci) and il{k2) (see Fig. The energy con- 
servation implies fl(ki) + Q{k2) = h(uji — cuf) = huj. 
The transferred frequency u is usually order of magni- 
tude smaller than tOij. 

Quite generally, the matrix element Mn depends on 
the two frequencies uji and Uf, the two momenta, either 
of two photons or two magnons, and on the polarizations 
Gi and Gf of the incoming and outgoing light. For the two 
most frequently studied cases of Aig and Big scattering, 
the polarizations are e.^ = (x + y)/ \/2, e/ = (x + y) /\/2, 
for scattering, ande^ = (x+y)/\/2, e/ — (x— y)/\/2, 
for Big scattering. For non- resonant scattering, which we 
consider, the frequency dependence of Mji only affects 
the overall factor, and the interesting physics comes from 
the momentum dependence of M^. The two magnon 
momenta can be parameterized as ki = k + q/2, k2 = 
—k + q/2, where q is the transferred photon momentum, 
then Mr = Mii{k, q). At g = 0, Mr oc cos kx + cos ky in 
Aig geometry, and Mr cx cos k^ — cos ky in Big geometry. 

Without the final state interaction, the intensity of the 
absorption of light RQ{q,uj) is obtained from the Fermi 



golden rule, and is given by 
Roiq.uj) (X Y.Ml(k,q) 



X 5[nLj-{n{k + q/2) + n{-k + q/2)].{l) 

Below we label by R{q, to) the full intensity, which in- 
cludes the final state interaction. 

In Raman scattering experiments with visible light, 
typical frequencies of the incident light are LOij / (27r) ~ 
510^^ Hz, and the corresponding wave-vectors qij ^ 
1.510^ m~^. Typical momenta k for magnons near mag- 
netic Brillouin zone boundary, which mostly contribute 
to Raman scattering, are of the order Tr/d, where d 
10^ ^''m is the lattice constant. Accordingly, k ^ 10^" 
m^^, three orders of magnitude larger than in which 
case one can safely set q = in Eq. 4]. The ap- 
proximation R{q,Lo) « R{q — has been used in all 
previous calculations of the Raman intensity 5, 6, 7]. 

In recent years, however, a new resonant inelastic X- 
ray scattering technique (RIXT) has been developed, 
which allows one to probe the intensity R{q,uj) at a fi- 
nite g H, 0) ^1 ■ The frequencies of X rays are tuned to 
the energies of atomic transitions and are several orders 
of magnitude larger than for conventional Raman exper- 
iments: uii/{2iT) ~ 210^* Hz, The corresponding 
momenta are then qij ~ 610^° m~^. These momenta 
are comparable to the typical magnon momenta k. As a 
result, RIXT technique allows one to measure the light 
intensity R{q, uj) at a finite q. 

A RIXT study of the undoped and weakly doped 
La2-xSrxCu04 and undoped Nd2Cu04 has recently 
been carried out by Hill et al ,1C|. In both undoped com- 
pounds, the authors scanned various momenta q, and 
at g = go = (0,7r) observed a sharp peak in i?(go,w) 
at 500meV. On deviations from go, the peak moves to 
smaller frequencies and broadens. The doping depen- 
dence of the peak at go is similar to that of the two- 
magnon Raman peak ~ both broaden with doping and 
disappear near optimal doping. Based on this similar- 
ity, the authors of IQ] argued that the peak at go is 
likely a two-magnon feature at a finite g. This is cor- 
roborated by the fact that the density of states (DOS) 
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of two non- interacting magnons has a peak at (0, tt) at 
3.5 J ~ bOOmeV. 

The presence of the peak in the DOS, however, does 
not necessary imply the peak in Ro(qQ,u!) as the matrix 
element M/j may counterweight the DOS effect. Besides, 
magnons do interact in the final state, and this interac- 
tion may substantially affect the profile of R{q, to) com- 
pared to that of Ro^q^uj). 

In this communication, we analyze R{q, lo) in Aig and 
Big geometries in two-dimensional (2D) antiferromag- 
netic Mott insulators (e.g., parent compounds of high- 
temperature superconductors). Like we said above, at 
q = 0, the Raman response in the Big geometry contains 
the two-magnon peak at w « 2.7J. In Aig geometry 
the Raman intensity vanishes at g = because Loudon- 
Fleury and Heisenberg Hamiltonians commute |^ |^ . We 
consider how both features are modified at a finite q. We 
show that the position of the two-magnon peak in the 
Big geometry evolves with q, and at go the peak is lo- 
cated at 2.4J ~ 3A0meV. We show that the intensity of 
the Big peak actually increases at a finite q and is the 
largest at qo and symmetry-related points. We further 
show that in Aig geometry, R{q, to) is non-zero at a finite 
q, and the intensity is again the largest at qo. We found 
that for q = qo, the intensities R{q,uj) in Aig and Big 
channels are equivalent both with and without final state 
interaction. 

Our results cannot be directly applied to the exper- 
iments by Hill et al |0 as they found the peak at qo 
in the geometry when the incident polarization is along 
c— axis. To obtain the matrix element in such geometry 
one would need to consider hopping between 2D planes. 
At the same time, our Aig and Big results show that 
the profile of the peak in R(qo,uj) is predominantly de- 
termined by the interplay between the two-magnon DOS 
and magnon-magnon interaction, while the momentum 
dependence of the matrix element does not play a 
substantial role near go- We therefore expect that the 
(identical) profile of R{qo,uj) in Aig and Big geometry is 
similar to that for c— axis polarization of the light, and 
the intensity is just stronger in c— axis geometry. 

The derivation of the matrix element Mn for x— ray 
scattering which involves high-energy photons in general 
requires a careful consideration of core atomic transitions 
\wL Il3 | . The outcome of this consideration is some effec- 
tive direct coupling between a:— rays and two magnons. 
In the following, we adopt an approach borrowed from 
Raman studies and derive the Raman matrix element 
by the same method as in Ref. 5], which we extended 
to finite q. Namely, we depart from the 2D half-filled 
Hubbard model with strong on-site repulsion U and the 
nearest-neighbor hopping t, expand the hopping term in 
vector potential A, and restrict with the on-site interac- 
tion term Hint — ~(e/?ic) jgA_q between the vector 
potential and the fermionic current. We then introduce 
antiferromagnetic long-range order parameter A (« U/2 
at U >> t which we only consider), and split quasipar- 
ticles into conduction and valence bands, with energies 



±-\/A2 + e^., where = —2t{coskx + cosky).'We then 
obtain from the model the interaction vertices between 
magnons and conduction and valence fermions. At large 
U/t, the dominant interactions are the ones in which 
fermion flips from the conduction to the valence band 
and vise versa after emitting a magnon. To keep calcula- 
tions under control, we extended the model to "large S" 
by introducing 25* fermionic flavors .12] . For large S, the 
corrections to our vertices are small, of order 1/S. 

Like we just said, for energies of the incident phonons 
well above Hubbard U, the matrix element Mfi comes 
from the core level atomic transitions rather then from 
the excitations across the Hubbard gap. Still, we believe 
that the derivation using "Raman technique" is useful as 
it shows the difference between the conventional Raman 
scattering with q = and a finite q scattering in dif- 
ferent scattering geometries. Besides, as we will see, for 
scattering near (0, tt) the profile of R{q,Lo) is determined 
predominantly by the structure of the DOS and magnon- 
magnon interaction, while the matrix element Mfi can be 
safely set to be a constant. This implies that the results 
below could be obtained with just a phenomenological 

Mr. 

The processes which contribute to the matrix element 
Mr in finite q Raman scattering are the ones in which in- 
cident light creates a particle-hole pair, which creates two 
magnons, either from two different lines in the particle- 
hole bubble, or from the same line, and then annihilates 
into outgoing light. The corresponding two diagrams are 
presented in Fig.^ [The diagram b has a "parasitic" con- 
tribution which for g = is canceled out by the additional 
diagram with an extra fermion- fermion interaction j5|. 
We verified that the same cancellation occurs also at a 
finite q]. Evaluating the two diagrams and adding them 
up, we obtain the matrix element in various geometries 
in the form 



MR{k, q) = -8i2 



2A 



4A2 



{[e„e}^ cos((7^/2) + e^ye) cos{qy/2)]{Xk+g/2fJ.-k+q/2 



+X-k+q/2l^k+q/2) " [^ixe}^ cos(fca;) -l- e^j^ej^ cos(fcj^)]x 



^ (Mfe+9/2/^-fc+g/2 + ^k+q/2^-k+q/2)}, 



(2) 



where 



1 

V2 



1, Xk 



and 7fc = (cos(fc^) -I- cos(fcj^)) /2. The magnon energy flk 
is given by ilk = Emax\/l - jI, where Emax = ^JS{= 
2 J), and J = t^/AU. At g = 0, the expression in Eq. ^ 
coincides with that in earlier works |^ . In particular, 
one can easily verify that MR^k, 0) vanishes in Aig geom- 
etry. At g 7^ 0, the matrix element is non-zero in all scat- 
tering geometries, including Aig. Furthermore, a simple 
trigonometric analysis shows that at q = qo = (0,7r), 



FIG. 1: The diagrams for the Raman matrix element 
Mii{k, q) Sohd and dashed lines represent conduction and va- 
lence electrons. The dash-dotted lines represent incoming and 
outgoing photons, and the solid wavy lines represent magnons. 
The full set includes an extra diagram (not shown) which can- 
cels out a parasitic contribution from the diagram b, see [3 • 




the values of Affl(fc,(7o) in Big and Aig geometries are 
identical. 

Without final state interaction, the scattering inten- 
sity Ro{q,ui) is given by Eq. In the top panels of 
Fig. 1301 we present the results for Ro{q,uj) for Aig and 
Big scattering. We see that in Aig geometry, the two- 
magnon peak gradually emerges as q increases, and pro- 
gressively shifts to larger frequencies. In Big geometry, 
the largest intensity without final state interaction is still 
at g = 0. At a finite g, the peak position shifts to lower 
frequencies. In Fig.^ top panel, we present the profile of 
Ro{qo,uj) (a dashed line) which, we remind, is identical 
in Aig and Big geometries. We see that Ro{qo,uj) has a 
sharp peak at ui — y/SEmax followed by the shoulder at 
u; = 2Emax , which is the maximum value of the energy of 
two magnons. The position of the peak at ^/SEmax can 
be traced back to the behavior of 5[uj — {^{k + qo/2) + 
Q{k — qo/2))], as -I- qo/2) + Q{k — qo/2) has an ex- 
tended van-Hove singularity at the set of fc— points for 
which Q{ko + (7o/2) + Sl(fco ~ Qo/'^) — V^Emax (see the 
bottom panel of Fig. The Raman matrix element 
M(k,qo) is regular and we found that its fc— dependence 
only weakly affects the form of Ro{qo,(^)- 

The Fermi golden-rule results for Ro{q,uj) are of lim- 
ited use, however, as two excited magnons interact in the 
final state. The magnon-magnon interaction is not small 
for S = 1/2, and in general substantially affects the scat- 
ting profile. We obtained the vertex for magnon-magnon 
interaction by a standard technique: we used the fact 
that the half-filled Hubbard model at large U reduces to 
the Heisenbcrg 2D antiferromagnet with nearest neighbor 
exchange, H = JJ2<i j> ' re-expressed spin opera- 
tors in terms of Holstein-Primakoff bosons a and f3 which 
describe excitations in the two sublattices, and restricted 
with only the interaction term a^/3^a/3 which describes 
multiple scattering of two excited magnons. This proce- 
dure is similar to the one used in previous works (see, 
G-g-, the new element is that now the total mo- 

mentum of the two magnons is non-zero. The relevant 
interaction part of the Heisenberg Hamiltonian is 

^' = Y. 'Y1 '^lik, 04+9/2'3-fc+?/2/5-i+«/2ai+g/2, 

k I 




FIG. 2: The scattering intensity in the Aig geometry along 
r - X - M - r direction, where F = (0,0), X = (vr, tt), and 
M = (0,7r). Top panel, without final state magnon-magnon 
interaction. Bottom panel - with the final state interaction. 
The largest intensity for the full R{q, cj) is at go = (0,7r). 



where the primes indicate that the summation is over the 
magnetic Brillouin zone, and 

r,ik,i) = -^x{ 

Ik^l {^J■k+q/2^J■l+q/2^J■-k+q/2^J'-l+q/2 + (^i ^ A)) 
+ lq {^J'k+q/2^l+q/2^-k+q/2^J'-l+q/2 + (A* ^ ^)) 

- 2 [ 7-fc+9/2 {f^k+q/2fk+q/2f^-k+q/2^-l+q/2 + (M ^ ^)) 
+ l-l+q/2 {^^k+q/2^J■l+q/2^-k+q/2^J■-l+q/2 + (P A)) 
+ lk+q/2 {tJ'k+q/2^l+q/2tJ'-k+q/2fJ'-l+q/2 + ^ A)) 
+ ll+q/2 (Afe+g/2/ii+g/2M-fc+g/2A*-i+9/2 + (A* ^ ^))] } ■ 

We computed the renormalization due to 
magnon-magnon interaction in the RPA ap- 
proximation, by projecting the interaction onto 
Mji{k,q) for different geometries, i.e., by approx- 
imating Tq{k,l) by BqMB,{k,q)Mii{l,q), where 

Bq = E[j^,{kJ)MRik,q)MR{l,q)/{j:[MUk,q))'. In 
this approximation, the total scattering intensity 



R{q, Lu) cx Im 



Ro{q,i^) 



l-BRo{q,Lj)' 



(3) 
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in the study of phonon-assisted optical absorption. In 
their case, a finite momentum q of two magnons is equiv- 
alent to a phonon momentum. Their effective matrix 
element for two-magnon vertex is different from ours and 
belongs to two-dimensional representation of the tetrago- 
nal group -D4/1. Still, our results are very similar to theirs 
- they also found that magnon-magnon interaction gives 
rise to a sharp peak in the absorption near 1.2Emax- The 
only distinction between our results and their is that in 
their case the peak moves to higher frequencies on mov- 
ing from qo towards (0, 0), while we found that Aig peak 
moves to smaller frequencies, and Big peak position re- 
mains almost unchanged 




FIG. 3: Same as in Fig|5| but for Big geometry. The intensity 
for the full R(q, uj) is again the largest for go = (0, tt). 



where 



Ro{q,uj) = ^ 



|M^(fc,g)|2 



k+q/2 



k+q/2) 



We set S* = 1/2 at the end of calculations. For q near 
qo — (0,7r), we obtained B w —3, i.c, the final state 
interaction is quite important. 

The results for the scattering intensity R{q, lj) renor- 
malized by the final state magnon-magnon interaction are 
presented in the bottom panels of Fig. |2 We see that 
both Aig and Big scattering are substantially sharpened 
up around q = qo, where the intensity is still at maximum 
in both geometries. 

We verified that the final state interaction does not 
break the equivalence between the scattering intensities 
in Aig and Big geometries right at go- In Fig- ^ top 
panel, we plot the full intensity R{qo, lu) versus frequency 
(a solid line). Comparing this result with i?o('?,w) (a 
dashed line on the same figure), we see that after the 
final state interaction, the peak gets sharper and more 
symmetric, and also shifts to a lower frequency, about 
1.2Ejnax- That the fully renormalized intensity R{q,Lo) 
has a sharp maximum at go is the central result of our 
paper. 

The final state magnon-magnon interaction at a finite 
q was earlier analyzed by Lorenzana and Sawatzky [l^ 




0.5 0.75 1 1.25 1.51 




FIG. 4: Top panel: the scattering intensities at go = (0, tt) 
with and without final state interaction (solid and dashed 
lines, respectively). Note that due to a final state interaction, 
the two-magnon peak gets sharper, becomes more symmetric 
and shifts to a lower frequency. Bottom panel: 2D constant 
energy contours for the total energy of two magnons -|- 
qo/2) + Q.{—k + qo/2). The extended van-Hove singularity is 
at fl{k + qo/2) +fl{-k + qo/2) =2^3.7. 

Our results for the full R{q, lo) partly agree and partly 
disagree with the experimental data for the inelastic reso- 
nance X— ray scattering of antiferromagnetic parent com- 
pounds La2CuOi and Nd2CuOi 10]. On one hand, the 
position and the profile of the peak fully agrees with our 
R{q,uj). Namely, the largest peak is at g = goj and it 
moves to smaller frequency and broadens on deviations 
from o^. On the other hand, we found, in agreement 
with [131 that due to final state interaction, the peak 
in R(qo,uj) shifts down to 1.2Emax = 2AJ ~ SAOmeV. 
This frequency is smaller than the experimental peak 
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frequency of 0.5eV « 3.6J, which is close to the peak 
position at qq in the absence of the final state interac- 
tion. The peak shifts to a higher frequency if we include 
quantum corrections to the magnon dispersion '17*1, ne- 
glected in our analysis. Another possible explanation of 
the discrepancy may be that the RPA-type analysis of 
magnon-magnon interaction overestimates the effect of 
the final state interaction. For c— axis polarization of 
light, as in the final state interaction involves the 

exchange between layers and generally should lead to a 
smaller down-shift of the peak position at (0, tt) from that 
fore non-interacting magnons |l8| . On the other hand, we 
didn't find the sharpest peak at (0, tt) for non-interacting 
magnons, e.g., in Big geometry, the peak at (0,0) is the 
sharpest. Magnon-magnon interaction substantially re- 
duces this last peak, and the one at (0, tt) becomes the 
largest. However, this only happens if magnon-magnon 



interaction is strong enough. 

To summarize, we analyzed the profile of the X-ray 
scattering intensity R{q, uj) at finite g in a 2D Heisenberg 
antifcrromagnet. Previous Raman scattering studies at 
g = identified the two-magnon peak in Big geometry. 
We found that the Big peak slightly disperses downwards 
at a finite q, and its intensity increases, reaching its max- 
imum aX q — qQ — (0, tt) and symmetry related points. 
Simultaneously, at a finite g, the intensity in Aig geom- 
etry also becomes non-zero, and the Aig Raman profile 
displays a two-magnon peak which disperses to larger 
frequencies with increasing q and reaches its highest in- 
tensity at go- We found that the profile of R{qo,uj) is 
equivalent in Aig and Big geometries. 

We thank G. Blumberg for useful conversations, con- 
structive criticism, and for sharing Ref. |lOj | with us. The 
research was supported by NSF DMR 0240238. 
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